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A mesoscopic evolution equation for an ensemble of mesoparticles follows 
after the elimination of internal degrees of freedom. If the system is composed 
of a hierarchy of scales, the reduction procedure could be worked repeatedly 
and the characterization of this iterating method is carried-out. Namely, a 
prescription describing a discrete hierarchy of master equations for the den- 
sity operator is obtained. Decoherence follows from the irreversible coupling 
of the system, defined by mesoscopic variables, to internal degree of freedom. 
We discuss briefly the existence of systems with the same dynamics laws at 
different scales. We made an explicit calculation for an ensemble of parti- 
cles with internal harmonic interaction in an external anharmonic field. New 
conditions related to the semiclassical limit for mesoscopic systems (Wigner- 
function) are conjectured. 
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I.- Introduction : reduction technique 



The study of models allowing a unified description of microscopic and 
macroscopic physical systems has a long history. The problem is related to 
microscopic superposition of states and its non-occurrence at macroscopic 
scales. Interesting responses and different proposition can be found in ref- 
erences [1-18] which are related to coherence destruction by different ap- 
proaches to the macroscopic level. Many of these theories are related to the 
original ideas developed by Landau [1] respect to the high density of states for 
macroscopic object, and high sensibility to external perturbation. Thus, any 
small perturbation produces an undefined (macroscopic) state and then mix- 
ture (i.e. decoherence) . Nevertheless, it is naive to think that the physics 
between elementary particles and macroscopic objects, for instance like to 
macro-molecules, can be described only for one mesoscopic theory. This is 
the case of the DNA-macromolecule which, in a first level, is composed of 
interacting atoms and finally becomes responsible for transmission of genetic 
information in a biological level. Other examples are some biological com- 
posites like to the hierarchical organization of tendon, bone, mollusk shell, 
synthetic composites [19,20] and others. 



The scope of this paper is the study of an iterating systematic procedure 
at different discrete scales of perception. The idea is simple, we start with 
an ensemble of elementary particles forming clusters around its mass-centers, 
then we eliminate the internal degrees of freedom. In this way the dynamics 
law for this reduced system is obtained and the procedure could be worked 
again to the next level. 

The general assumptions to construct the hierarchy of master equations 

are: 
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(i) The system is composed of a hierarchy of scales, or levels, 
where we can recognize different architectures of clusters. 



(ii) There are intercorrelation between levels. Specifically, a given 
set of clusters make-up a cluster in the next scale. 



(iii) The dynamical requirement (equations of motion), in a given 
scale, is depending onto the above one. 



Assumption (i) and (ii) are related to geometrical aspects, and (iii) pos- 
sibilities the obtention of dynamical laws from more fundamentals scales. In 
fact, (iii) is related to the usual belief that phenomeno logical laws can be 
explained from fundamental models. 



The paper is organized as follows. In this first section, we discuss briefly 
the reduction method applied to an ensemble of generic systems. In section II 
we deal with the iterating method and finding the master equation describing 
the dynamics, at different discrete scales, for interacting mesoparticles. In III 
some examples are briefly studied and we investigate the question related to 
the existence of invariant-systems under reduction procedure. In section IV, 
we use the Wigner function to explore the semiclassical limit for the quantum 
evolution equation. The mesoscopic term, related to the internal degree of 
freedom, requires new conditions aside the usual one related to the optical 
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geometric limit. Conclusions and discussions are presented in the ending 
section. 

Now we reviewed briefly the reduction procedure [21-27] which will be 
used in the next section. Consider the interaction between a system S and 
other R, with many degrees of freedom, and the evolution equation for the 
complete system 

d t p = Cp. (1) 
We assume that the Liouville-von Neumann operator is decomposed like 

£ = C s + Cr + Ci (2) 

where Cj denotes the interaction term. Consider the projector operator P 
acting onto the total density operator p (or the space density distribution in 
the classical case) : Pp = p e R Tr R (p) where p e R denotes the equilibrium state 
of the system R. In the classical case, the partial-trace-operation is replaced 
by an integral over the phase-space of R. As usually, the projection operator 
P satisfies [22] : 

PC R = C R P = 0, PC S = C S P, P£/P = 0. (3) 

In this way from the evolution equation (1), for the complete system, and 
projecting on the space spanned by P and Q = 1 — P, one obtains the 
equation for S : 

d tPs = C sPs + Tr R C T f dre CT C lPR p s {t - r). (4) 

Jo- 

To obtain the above equation the usual initial condition p(0) = p e R ps was 
assumed. Equation (4) is exact and cumbersome because the integral term is 
dependent on the history of S. Expanding to second order in the interaction 
term one obtains the equation 

dtPs = C sP s + f dTid&WnPsit), (5) 

J 
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where (o}^ denotes the partial trace operation Tr R (op^), and the symbol £ o' 
means an element of the space of operators. Assuming the usual memory 
loss property or Markov approximation : 



with 7 a positive parameter and Cis an operator acting on ps, then we obtain 
from (5) and (6), the evolution equation for S 



The explicit verification of the properties : ps = p$, Trps = 1 and ps > 
(positivity) must be carried-out always. Remark that currently an equation 
like (7) is related to decoherence. Specifically, the reservoir R changes any 
system pure states to mixed. An approximated case, where Cjs = £s, with 
decoherence and without dissipation can be found in [8,9]. For a criticism to 
the reduction procedure see for instance [10,12] where a completely integrable 
system was considered. Dissipating effects are treated in [4] (non-linear equa- 
tion) where other techniques were considered. 



II.- Reduction procedure: hierarchiy of mas- 
ter equations 



In this section we use the reduction procedure sketched in section I, 
including a coordinate change to the center of mass, and we obtain the 
evolution operator at the next scale. Namely, we formalize the procedure 

£(n) _^ £(n+l)_ go at gcale 

n, we start with interacting particles in an 
external field. This set contains an architecture of jV^ -1-1 ) clusters which are 



(6) 



d t p s = £>sPs + itfsPs- 



(7) 
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recognized by using physical constraints (assumption (i)). For instance they 
could be N^ n+1 ^ molecules in an electric field or interacting macromolecules. 
After the cluster recognition, we consider a coordinate change to the center of 
mass (assumption (ii)), eliminating the internal degrees of freedom for every 
cluster (tracing-out-technique). 

Formally at scale n, the equation of motion for the density distribution 
in the classical case, or the density operator in the quantum case, is 

dtf) {n) = £(n) p (n) ; (g) 

where denotes a linear operator constructed, for instance, by elemen- 
tary Liouvillian (or von Neumann) operators like £/o = {/, o}, (or CfO = 
(l/ih)[f, o]). Where the symbol{o, o} ( or [o, o]) stands for the usual Poisson 
brackets (commutator). 

In (8), the index n becomes related to the discrete scale, and the idea 
is to obtain the evolution equation, at next scale, by constructing the new 
operator £( n+1 ) from the dynamics at scale n (assumption (iii)). In this way, 
the reduction procedure has technically two steps: 



(a) A coordinate change to the center-of-mass of every cluster. 

(b) Elimination of internal degrees of freedom by assuming 
loss- memory-effects (i.e. internal complexity). 



Therefore, the reduction — * jV^ 1 ) gives an equation similar to (8), 
where the Liouville-von Neumann operator is determined using steps (a) and 
(b). To find the new evolution operator, we assume the decomposition : 

£ (n) = £ (n) + £ M (Q) 
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Namely, a kinetic part depending on momentum and other depending on 
position. 



To begin, we explicitly consider the first reduction procedure because it 
contains all the basic ingredients for further iterations. Namely, we consider 
the reduction — > where the index n = stands for an elementary 
set of interacting particles forming clusters. 

Let Xj( a ) be the position of the particle j (integer) in the cluster a (integer) 
where 1 < a < N^. Consider the transformation to the center of mass y a , 
of the cluster a, given by 



where rj( a ) denotes the relative distance with respect y a . Consider the inter- 
acting internal potential Ut, 



where q is a coupling parameter and the summation rules on all indices and 
no self-interactions, or repeated indices, are assumed. Moreover, like j(a), 
the term k((3) denotes the particle k in the cluster at yp. From (10), the 
potential transforms like 



U t = E ( U (y° - VP) + (O(a) - r m )U'(y a - y p )) + F(r) (12) 



where a first order multipolar expansion, in the internal coordinates, was 
assumed and the symbol U' denotes the first derivative. Moreover in (12), 
F(r) denotes the linear terms depending only on the relative coordinate and 
related to internal interaction in the cluster. Now we define the coupling 
parameter, between clusters, Q 1 ^ and the moment d a> p like 



(10) 



(11) 



Qi = H 9j(a) ' = E 1j(a)( r j(a) ~ ^■( / 3))- 



(13) 




7 



and from equation (12) and (13) the internal interacting potential becomes 
Ut=E QiU{y a - y p ) + dafU'iy* - y p ) + F(r). (14) 

a>(5 

Similarly, for an external field acting on every particle and given by 

V T = ^2qj{a)V(x j{a )) (15) 
<x,3 

and denning the new coupling parameters and dipolar distribution by 

Qoc = m ° = Qj(a) r j{a), (16) 

J (a) 

then the potential (15) can be written, at first order in the relative coordi- 
nates, as 

V T = Y,Q a V(y a )+m a V'(y a ) (17) 

a 

and like (12) the symbol V denotes the first derivative. Remark that no term 
like F(r) appears in this case. From (14,17), and the transformation for the 
kinetic term which is form invariant, the complete Hamiltonian becomes 

H = H s + H R + Hj (18) 
where the explicit form for the Hamiltonians are 

2 

H * = Y,£r + Q*v(y a ) + E - vj) (w) 
Hi = J2 V(y a )m a + J2 u '(y« - vp)**,? ( 20 ) 

a a>/3 

and Hr is the contribution due only to internal coordinates (r, f). In equation 
(19) the term [i a denotes the total mass of the cluster a. Recall that it was 
always assumed a first multipolar order expansion in the internal coordinate. 
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In this way, the Hamiltonian (18) has the structure worked in section I. 
Since the internal variables d a $ and m a are assumed independent and they 
have the loss-memory property then , from section I, the evolution operator 
for the ensemble of mesoparticles with position y a is 

= £<°>(y,p) + £ 7 « (£$(y a )) 2 + £ 7 « (4g(y a - y,)) 2 (21) 

a a>/3 

where we have for every component the expressions : 

£(°\y,p)o = {H s ,o} (22) 

(£§(y a )) 2 o = {V'(y a ),{V'(y a ),o}} (23) 

(4g(y a - ^)) 2 o = {tf'(y a - y/J ), {[/'(y Q - y/J ), o}} (24) 

namely a double Poisson brackets, or double commutator in the quantum 
case. 

In (21) the parameters 7 are related to white- noise type correlations be- 
tween internal variables (dipole moments) for every mesoparticle . This is 
the Markovian approximation where memory effects are ignored. Expression 
(21) gives us the evolution operator £^ for mesoparticles of coordinates 
(y a ,Pa), where internal degrees of freedom were eliminated. We noticed that 
assumption (iii), of section I, is in accord with our deduction because the 
evolution operator at scale n = 1 was deducted from this one at scale n — 0. 
Similar equations for one, or two mesoparticles, were also discussed in [28- 
30]. The idea to use internal degrees of freedom as an internal environment 
are also discussed in [31]. 
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The above evolution operator (21) is related to the first elimination of 
internal degrees. Nevertheless, if we can recognize a second structure of clus- 
ters, we can eliminate new internal degrees. Evidently, this will be possible 
only if the Markovian approximation is valid. Let C^(xj,qj) be the evo- 
lution operator at scale n, which includes kinetics and potential terms like 
(9), and consider a cluster recognition with center of masses at y a where 
1 < a < N( n+1 \ Then the formal first order multipolar expansion with 
respect to the relative coordinates (10) is 

C <»> = £<»> (y, Q) + £ dCi y> Q) d a + F(r) (25) 

where the formal derivative stands for the first order expansion, and d a are 
linear function of the internal degrees. Moreover, Q denotes some re-defined 
coupling parameters. Then still we have a situation similar to this of section 
(I), and if we assume the loss memory effect then the new evolution operator 
for the TV*™" 1 " 1 ) mesoparticles will be 

£ (n + D =£W (y ,Q)+£ 7 W ( g£ ^ g) ) 2 . (26) 

At this point some remarks related to the iterating equation (26) are: (a) 
The formal derivative stands for the first order multipolar expansion around 
the center of mass of every mesoparticle, or cluster, (b) Expression (26) is 
valid in the quantum or classical case, where C becomes related to a set of 
elementary commutators or Poisson brackets, (c) Space structure (lattice, 
fractal, disordered, etc.) is contained in C°. It decides the criterion for the 
cluster recognition, (d) Remark that the free particles case (£y = Cu — 0) 
is a trivial form-invariant-example under reduction procedure. 
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III.- Examples 



We will now examine briefly some examples related to the reduction pro- 
cedure discussed in section II. Explicitly, we consider an ensemble of harmonic 
interacting particles with anharmonic external fields. In fact, the structure 
of the evolution operator becomes invariant after some reductions (aside of 
some redefined coupling parameters). The search for such an invariance was 
also investigated in [3] for a two-parameter model of decoherence, where the 
equation of motion for the center of mass, is formally identical to the equation 
for the microscopic constituents. 

Consider an ensemble of particles with harmonic interaction, in a nonlin- 
ear external field. The internal and external interaction operators are given 



where Kij are positive constants and A a coupling parameter. The expansion 
around center of mass gives 



by 




(27) 




(28) 
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E( 1 /2){^(^ W -^(/3)) 2 ? °}+ 



+ 2 S{^i(a),fcOJ)(ya - yp)(rj(a) - r m ), o} + F(r) 



(29) 
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and for the external component 

4 0) ° = a' E {vl °} + 3A E {vfow. °} (30) 

where if and A' are re-defined coupling parameters. The reduction procedure 
(section II) leads to the evolution operator 

£ (D = E{|^ + A'j£, o} + E - ^) 2 , o} + 

+ E <p {(v« ~ vp) diva- y/») ,°}} + E 7^ {y«, R °}} (3i) 

which gives us the evolution equation for the set of mesoparticles. 

At this point we have an interesting result, a second reduction process, 
makes invariant the (internal) evolution operator. In fact, the only changes 
are related to the redefinition of the coupling parameters and mass. The 
same is true for the external anharmonic term in (31), which becomes invari- 
ant after three reduction process. This seems a general fact related with the 
formal derivative in the expression (26) for the evolution operator. So, for 
systems interacting algebraically (i.e. V, U ~ x n ) the reduction procedure 
seems invariant after a number finite of steps. Namely, the laws of evolu- 
tion become the same at different scales of perceptions. Nevertheless, the 
reduction procedure must be stopped when there is not loss-memory-effect 
and then this process does not can be continuously carried-out. This is the 
case for instance for a system with a finite number of constituent. Also, we 
noticed that geometric aspects must be considered at every reduction, and 
some interesting candidates for a such invariance are elementary excitations 
in fractal structures. 



12 



IV.- Wigner function and the classsical limit 



Since (21) is also valid for quantum system making the appropriate changes, 
it is instructive to study the semiclassical limit using the Wigner function. 
In fact, we shall find that the semiclassical limit needs some new conditions 
because the decoherence term related to the parameter 7. 



The Wigner function p w defined by the Fourier transform of the density 
operator in coordinate representation is given by 



Pw (x, P,t) = lJ d V e^ h p(x - 77/2, x + 77/2, t). 



(32) 



Where p(x, y, t) is the statistical operator in coordinate representation. Then, 
from equation (21) and keeping by simplicity only the external potential term, 
the evolution for the Wigner function becomes 



d t p w = {H s ,p w } + J 



dV 

dx 



dV 



Pv 



+ (h 2 /24) 



'd 3 Vd 3 p h 



-2 7 



„ ^x 
d 2 V d 4 V d 4 p 



+ 



0{h 4 ). 



(33) 



dx 3 dp 3 dx 2 dx 4 dp 4 

The first two terms are the classical operators (22-23) and the other are 
related to the quantum contribution. So, the mesoscopic term related to 
the parameter 7 gives new quantum corrections. The usual semiclassical 
approximation, when 7 = 0, is the well-known relationship 



dVdp u 



dx dp 



»ft 2 



d 3 Vd 3 Pu 



dx 3 dp 3 



(34) 



and related roughly to the optical geometrical limit. It must be noted, how- 
ever, a similar condition related to the mesoscopic term from (33) this con- 
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dition is (7 7^ 0) 



'd 2 vYd 2 Pv 



dx 2 J dp 2 



d 2 V d 4 V d A Pn 



dx 2 dx A dp A 



(35) 



Namely, it explores even derivatives in the potential. 



Finally we note that the deduction of the evolution equation (33), for the 
Wigner function, was obtained assuming (integration by part) 

ihd x p(x -r],x + ^e^lJpS = 0. (36) 

Namely, the vanishing of the correlation term at infinite. This requirement 
is not always verified, for instance, consider the states p = ^{y)*^{x) where 
the wave function is i/j — sinkx which does not satisfied (36). Nevertheless, 
the contribution due to the decoherence in (21), produces a fast annulment 
of the off-diagonal terms in the statistical operator for short-range external 
potential. In this way, the condition (36) can be satisfied for times greater 
that the decoherence time, and given a solid support to the evolution equation 
(33) for the Wigner function. 



V.- Conclusions and discussions 



We have considered a hierarchy of master equation describing the evolu- 
tion, at different scales of perception, for ensembles of mesoparticles. Explic- 
itly, the equation (26) gives us formally the evolution operator £( n+1 ) from 
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this one at scale n (assumption (iii)). Its deduction requires a systematic co- 
ordinate change to the centers of mass, defined by some physical constraints, 
and the elimination of internal degrees of freedom is carried-out assuming 
loss-memory-effects. This Markovian approximation is not always valid and 
then, in such a case, the reduction process must be stopped. The first re- 
duction procedure was carried-out explicitly for an ensemble of elementary 
components (21) 

It must be noted that assumptions (i— iii) , of section I, are the basis where 
our reduction procedure was developed. Namely, they possibility the obten- 
tion of a hierarchy of master equations at different scales of perception for 
cluster of mesoparticles. Some similarities between (i— iii) and those used in 
the architecture of complex synthetic assembles would be found in reference 
[19,20]. 



On the other hand, some important points related to the deduction of 
the evolution equation (21) and (26) are : 



(a) The first multipolar order expansion, in the interaction term, tell us 
that the asymptotic limit t — > +oo must be carried-out carefully at different 
scales [29]. 

(b) Decoherence effects at every discrete scale appears usually related 
to the reduction technique (section I). Thus, decoherence at macroscopic 
level is due to the internal complexity of every macroparticle. And quantum 
superposition is turned into statistical mixture . 

(c) The search for invariant systems was carried-out explicitly for a model 
composed of interacting oscillators with an anharmonic term. It seems that 
other invariant systems could be found. 



15 



(d) It was assumed that the internal moments, like to d or m in (13) and 
(16), are random and independent. Obviously this is not easy to prove and 
we have only assume that behavior. The statistical independence, between 
these random variables, is a simplification related with our calculations. 

(e) The deduction of a general equation like (7) was carried-out assuming 
special initial conditions (p(t = 0) = p e R ps)- In our specific case of section II, 
these conditions not necessarily hold. More explicitly, the internal interaction 
cannot be switched-up arbitrary. It is an open problem to prove the validity 
of our procedure in this case. 



To explore possible applications of our prescription, we can consider 
trends like wavematter, currently studied theoretically as well experimen- 
tally. After all, interacting atoms in external fields could be considered as 
mesoparticles. Specifically, with laser cooling techniques it becomes possible 
to cool atoms so that the quantum nature of atomic center of mass motion 
becomes important [32,33]. Also, it can be interesting note the growing inter- 
est in new mechanisms to break Anderson localization in disordered systems 
[34]. Particularly, there is the controversy about the possible enhancement 
of the localization length for interacting particles (TIP) [35]. This suggests 
considering the behavior of mesoparticles in disordered systems. Namely, an 
equation like to (21) with external random potential. After all, localization is 
a phenomenon related to coherence which is loss due to internal complexity 
for mesoparticles. A more detailed treatment of these points, and further 
physical applications will be given elsewhere. 
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